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Abstract The quantum evolution equation of Loop Quantum Cosmology (LQC) - the quantum 
Hamiltonian constraint - is a difference equation. We relate the LQC constraint equation in vacuum 
Bianchi I separable (locally rotationally symmetric) models with an integrable differential-difference 
nonlinear Schrodinger type equation, which in turn is known to be associated with integrable, discrete 
Heisenberg spin chain models in condensed matter physics. We illustrate the similarity between both 
systems with a simple constraint in the linear regime. 
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^vq , 1 Introduction 

VO, 

Loop Quantum Gravity (LQG) [1] aims to describe the physics of spacetime at Planck- 
length scales. In the last few years, research in this area has progressed towards predic- 
tions that area and volume observables are fundamentally quantized. In addition, LQG 
has been explored in symmetry reduced models, which take in consideration simple 
cosmological spacetimes, a research area known as Loop Quantum Cosmology (LQC) 

[2]- 

The quantum evolution equation of LQG - the quantum Hamiltonian constraint - is 
a difference equation for the wavefunction (i.e., a discrete recursion relation), instead of 
a differential equation like the Wheeler-de Witt equation (see, e.g., Ref. [3]). A series of 
consistency studies indicates that solutions to the quantum evolution equation, under 
certain conditions, do not lead to singularities: a minimum observable volume leads 
to a "bounce" near classical singularities and quantum states evolve deterministically 
through those singularities (e.g., [4], [3]). Recently, the nature of the bounce and its 
robustness have been addressed under several approaches: numerically; by use of ap- 
proximation methods ("effective equation" techniques); or by the analysis of simplified 
exact analytical solutions (for reviews, see, e.g., [6], [5] and [7] and references therein). 

However, as symmetry-reduced versions of LQG, LQC models assume certain results 
directly from full LQG. For instance, the smallest nonzero area eigenvalue of LQG is 
the assumed step size (the so called "area gap" ) of the LQC difference equation. It has 
been shown that the Wheeler-de Witt equation can be made to agree with LQC (in the 
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case of a k = 0, A = FRW cosmological model, coupled to a massless scalar field), 
as the area gap diminishes. However, the approximation is not uniform in the chosen 
interval of "internal time" , and in fact if the area gap is set to zero, LQC does not admit 
any limit, being an intrinsically discrete theory [4]. An elucidation on the nature of the 
solutions of the discrete Hamiltonian constraint in more general cases is necessary. 

Here, we propose a technique which may be able to address these questions more 
generally, by observing that solutions to the fundamental equation of micromagnetism 
in condensed matter have been explored in the context of similar difference equations. 

2 LQC and Nonlinear Schrodinger equations 

Generally, the understanding of the magnetization vector field dynamics of a submicron 
magnetic particle is made at at the continuum approximation, according to the Landau- 
Lifshitz equation (LL; also known as the Landau-Lifshitz-Gilbert equation with the 
inclusion of damping) [9]. It describes the local precessional motion of the magnetization 
vector under an effective magnetic field which represents several distinct interactions 
amongst the spins. Being a highly nonlinear vector partial differential equation, the 
LL equation is generally solved by using numerical techniques (analytical solutions 
can only be obtained in very few cases; see, e.g., [10]). In particular, the dynamics of 
classical Heisenberg ferromagnetic spin chain models is governed by the LL equation. 
Such models are interesting nonlinear systems, displaying a rich diversity of dynamical 
properties, including soliton spin excitations [11]. 

It has been realized, however, an interesting relation between such models and a fam- 
ily of completely integrable nonlinear Schrodinger (NLS) equations [12]. For instance, 
it can be shown that, given the one- dimensional isotropic ferromagnetic spin chain clas- 
sical Hamiltonian H = — J2n S n • S n+ i (wherein the resulting LL equation of motion is 
dS/dt = S x d 2 S/dx 2 ), there is an associated cubic NLS equation, which is completely 
integrable, given by i(dq/dt) + (d 2 q/dx 2 ) + 2|g| 2 g = 0. These results were obtained 
by the use of integrability methods such as the gauge equivalence and the geometric 
equivalence (space curve) formalisms ([13] and references therein). However, a more 
realistic description of ferromagnetic systems is a formulation of the dynamics at the 
lattice level, instead of the continuum limit. The integrability of such discrete systems 
has also been addressed under the methods above, suitably adapted for discretization 



In the present paper, we present a simple example which makes use of the above 
mentioned results of discrete Heisenberg ferromagnetic spin chain models in the con- 
text of difference equations of LQC models. We will focus on a simple LQC evolution 
equation in vacuum Bianchi I separable, LRS (locally rotationally symmetric) models, 
which are the simplest anisotropic model (with only two degrees of freedom), particu- 
larly we shall focus only on the resulting one-parameter, discrete recursion relation that 
arises from this model [15]: 



where q(m) (m, real numbers) is related to the coefficients of the wavefunction t m ^ n , 
which in turn are rescaled versions of the original wavefunction, normalized by a factor 



[19], [14]. 
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of the world volume of each basis state [15]. A is a separation parameter that arises in 
the scheme, and the function f(m) (f(m) = \ jm for m > 1 and f(m) = form m — 0) 
is a simplification of the function d(m) in the last Ref. of [15]. 

In parallel, we introduce the following completely integrable differential- difference 
nonlinear Schrodinger (NLS) equation (see, e.g., [14]): 



.dq(n) 
i- 



<jf (l + \q(n)\ 2 ) [q{n + 1) + q(n - 1)] - 2q(n), (2) 

where n represents an integer sequence of points. We parametrize the factor (1 + |<?(n)| 2 ) 
in that equation by the operator a, and write the NLS equation as: 

Qq(n) — q(n + 1) + q(n - 1), (3) 

where the operator Q is given by: 



Q = \ 

a 



■ d n 
*dt +2 



(4) 



We also rewrite the LQC equation, Eq. (1), as: 

fiq x (m) = q x (m + l)-q x (m-l), (5) 

with T = A/ (to). Whether the LQC and NLS discrete equations (Eqs. 3 and 5) can be 
made mathematically equivalent in form is a relatively straightforward question, as far 
as one assumes some restrictions, whose validity on physical grounds is, nevertheless, 
unclear at this point. We shall address possible implications as we proceed. Eqs. 3 and 
5 formally differ from each other in the following points: (i) the nature of operators Q 
and J 7 , (ii) the nature of indexes to and n (iii) the existence of a family of equivalent 
equations parametrized further by the A in the LQC equation, and (iv) the existence 
of a plus sign instead of minus sign in the right hand sides of these equations. Let us 
analyse each of these items separately. 



3 Linear sector 

First, equivalence between operators Q and T impose the restriction that the wavefunc- 
tions q\(m) should be of the form 

q x (m) = e WM-*o] j ( 6 ) 

where k x ~ aX defines a fiducial wavenumber; X(m) ~ f{m)i a fiducial sequence 
"length" , and $o = 2, a fixed phase. We notice here a similarity between that restriction 
and the choice of basis elements in Ref. [16], 

b k ( m ) = e ikexp (-&\ (7) 

where M a fiducial length. Such basis is adopted in that work in order to find pre- 
classical solutions to the separable LQC equation (Eq. 5). Pre-classicality is a con- 
cept introduced in order to discard unphysical solutions to the Hamiltonian constraint 
equations. These are solutions with highly oscillatory behavior arising far away from 
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the corresponding classical singularity; in other words, pre-classical solutions are those 
which the coefficients of the wave function at large volume should not oscillate when 
small changes to the volume (e.g., by amounts of Planck size) are made (for a critical 
review on pre-classicality, see Section 5.18 of the first Ref. of [2]). The idea in Ref. [16] 
is to write a generic solution to Eq. 5 as an expansion of the basis (7), which enables 
oscillatory behavior for small values of the index parameters, but at the same time 
singles out only smoother behavior for larger values of the parameters. We note that 
a choice of X(m) in Eq. 6 can be made to agree with the exponential factor in Eq. 7, 
generally, as: 

,w x <•/ x 1 / 2m. „ 2m . . 

X(m) ~ f(m) ~ — exp ( — — ) ~ 1 - — , (8) 
v ' v ' m v M ' M w 

which implies that M is now a function of m, 

2m 2 

M(m) ~ , m > 1, (9) 

1 — m 

in order that our restricted wavefunctions (6) allow pre-classical solutions in the form 
considered in Ref. [16]. It means that each m in the sequence is associated with a 
given fiducial M. Other arbitrary forms can evidently be assumed for X(m), which are 
completely independent on our assumptions concerning the equivalence between the 
LQC and NLS discrete equations (Eqs. 3 and 5), since for the moment such equivalence 
in terms of the operators Q and JF only restricts the form of the wavefunctions q\(m) 
in accordance with Eq. 6. 

Second, concerning the nature of the indexes, we only point out that pre-classical 
solutions to Eq. 5 are known for both when m is an integer or not, through the method 
of generating functions (c.f. [15]). For definiteness, we analyse only integer index param- 
eters, in accordance with the integer index parameters of the NLS equation formulation 
(Eq. 3). Third, concerning the parametrization of solutions according to A, we assume 
they enter our formulation through the restriction 6, which can be freely varied for each 
family of solutions dictated by A. Hence, we will assume that A is fixed a priori in order 
to proceed with the equivalence between Eqs. 3 and 5. 

Finally, we address freedom of signs (plus, minus) involved in the right hand sides 
of Eqs. 3 and 5. For that, we briefly outline the method of space curve formalism 
(mentioned in the previous section), which allows one to map continuum Heisenberg 
spin-chain models onto the NLS family of equations. Our analysis was based on the 
results obtained by the discrete mapping procedure of Daniel and Manivannan [14], 
motivated by the extension of the space curve formalism for discrete cases by Doliwa 
and Santini [19]. The reader is referred to these papers (and references therein) for 
details. 

The idea arises from well-known developments in mathematical physics, connecting 
the differential geometry of submanifolds and nonlinear partial differential equations. 
In essence, one looks for a fundamental geometric characterization of the motions of a 
curve in a submanifold which selects, among all possible dynamics, integrable dynamics 
[13]. The particular method for the problem at hand begins by considering a discrete 
curve on a sphere, represented by a sequence of points and a set of three orthonormal 
basis vectors associated with each point (defined by a position vector with respect to 
a reference point at the center of the sphere). One also specifies adequate transition 
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matrices from on basis to the next basis, in this case, these matrices can be combined 
into a rotation matrix depending only on two angles (9 n ,<f>n), namely, involving two 
of the basis vectors at subsequent points of the curve. The evolution of the basis is 
described by a matrix equation, and by defining a shift operator along the curve, as 
well as a compatibility condition for this operator, one reaches at a set of coupled 
evolution equation for the angles 9 n , (f> n . One then maps the Hamiltonian of the discrete 
chain model to be studied and its discrete equation of motion onto the discrete curve 
by a direct identification of the spin vector to one of the unit basis vectors (the radial 
vector). Such association leads to a completely equivalent set of coupled differential- 
difference equations for the evolution of the two angles 9 n , <p n . The question at this 
point is to verify whether such coupled equations are integrable by finding a suitable 
transformation to a known integrable differential-difference equation. 

By the method outline above, Daniel and Manivannan show that the completely 
integrable differential-difference NLS equation (Eq. 3) is associated with a classical 
equation of motion for the discrete Heisenberg spin chain, the Ishimori's model (based 
on Ishimori's partial differential equation, c.f. Ref. [18]), which is a classical discrete 
equation of motion, given by: 



^S n (t) — 2S n x 



S n +i S n _i 



(10) 



1 + S n • S n _|_i 1 + 

where S n = (S x , S y , S z ) are classical 3D spin vectors, with the constraint that |S n | 2 = 
should be kept as a constant. 

Now returning to the question of the freedom of signs involved in the right hand sides 
of Eqs. 3 and 5, by a detailed inspection of the equations of the mapping procedure 
outlined above, the LQC difference equation can also be mapped to Ishimori's model, 
inasmuch the coordinates are simply counter-rotated by a substitition sin <p n — > — sin <p n 
in the matrix given by Eq. (6b) in that paper, in order to account now for the minus 
sign in the difference equation (5). Therefore, the LQC difference equation presents a 
completely equivalent discrete spin chain equation counterpart, mathematically equiv- 
alent to the Ishimori's model, and therefore, under the conditions previously stated, it 
is integrable. 



4 Conclusions 



It is interesting to note that the equivalent counterpart of the complete Ishimori's 
equation [18] (spin-one field model in the plane) is given formally the Davey-Stewartson 
equation [20], which describes modulated nonlinear surface gravity waves, propagating 
over a horizontal sea bed (3D packets of mechanical surface waves) 1 . We speculate to 
what extent these models are somehow able to mimmic other classes of reduced LQC 
models, and if so, what they are telling us about the propagation through classical 
singularities. 

In summary, we have raised the possibility of relating the LQC constraint equation 
in vacuum Bianchi I separable models with an integrable differential- difference nonlin- 
ear Schrodinger type equation, which in turn is known to be associated with integrable, 



1 In 1+1 dimensions, this equation reduces to a simple nonlinear Schrodinger equation, Eq. 3 of the present 
paper. 
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discrete Heisenberg spin chain models. If this association proves to be correct at some 
level, then it is encouraging to have integrable equations at our disposal as LQC toy 
models, specially in the nonlinear regime. Here we have only pointed out an illustration 
in the linear regime. This proposal also opens the possibility of studying the propaga- 
tion of quantum gravity reduced wavefunctions models using spin chain representation 
models, which are amenable for computer simulations. Several already available numer- 
ical integrators of micromagnetism can be adapted to study these models. However, a 
deeper understanding of the physical meaning of "Ishimori-like quantum gravity mod- 
els" , specially the meaning of the vector spin, are still unclear at this point, and research 
in this area is under way. 

Acknowledgements The author acknowledges M. Bojowald for useful discussions, and the referees for rec- 
ommending some improvements in the exposition. 
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